Spin Splitting in Quantum Wires with Perturbed Axial Symmetry 
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We investigate theoretically perturbations to the confining potential capable of lifting spin de- 
generacy in axially symmetric quasi-one-dimensional electron gases with the spin-orbit interaction. 
The role of two diflFerent types of perturbations breaking axial symmetry of the system: of non- 
electromagnetic and electromagnetic origin is analyzed. We found that these two types of pertur- 
bations have a fundamental distinction in their effects on energy spectra. The influence of the scale 
of perturbation upon the value of spin splitting is investigated and it is shown that under certain 
conditions such splittings can be observed in the experiment. 
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I. INTRODUCTION 



Spin splitting in energy bands of quasi-one-dimensional 
and two-dimensional electron gases (QlDEG's and 
2DEG's, respectively) in zero magnajtiR-.£jG|ld has been 
studied (p±.pnsijHy-]3oth theoreticallylila tjQtJ and experi- 
mentallyEErEjild'Ej'Ej over the two past decades. Such an 
interest was partially produced by the prospect of cre- 
ating new electronic devices which exploit spin-polarized 
transportcJtJ. The appearance of spin splitting in energy 
spectra of realistic heterostructures, which are usually de- 
scribed by QlDEfG's and 2DEG's ipedels, was ascribed 
mainly to RashbaO and Dresselhausta terms in the corre- 
sponding Haniiltonians. Such Hamiltonians usually pos- 
sess two distinctive features: (i) as parameters varied 
in theoretical constructs they contain merely coefficients 
placed before splitting-inducing terms in order to achieve 
an agreement with the experimento and (ii) they include 
confining potentials for effective reduction of the dimen- 
sions of the system. While a thorough analysis of the 
effect of different splitting-inducing terms upon energy 
spectra has been carried out and their role understood, 
the choice of confining potentijals is still limited|«iainly to 
simple (and model) hard-walH and parabolicBLna'H ones. 
Such a choice is conditioned by the fact that the cor- 
responding non-perturbed problems are often explicitly 
solvable. At the same time this may give rather simpli- 
fied picture of real energy band structure of the model 
systems. The common step to avoid this shortcoming is 
to change the form of a potential slightly, for example, by 
accounting for non-parabolicity terms. Here, however, we 
do not touch upon this and related issues which are un- 
likely to seriously influence the general spectral picture, 
but introduce confining potentials which are themselves 
capable of giving rise to spin splitting beyond Rashba 
and Dresselhaus mechanisms. 

In this paper we investigate an infinite QIDEG pro- 
duced out of 3-dimcnsional system by means of a 2- 
dimcnsional confining potential. The initial confining 
potential is assumed to be axially symmetrical and have 
the electromagnetic nature. Axially symmetrical quan- 



tum wires are rather easy to investigate theoretically 
and, what is more important, they are not just a the- 
oretical mO|djel, but may be prepared practically (see, 
for instancolj, and references therein). We also neglect 
Coulomb interaction between electrons, which is quite re- 
alistic given that it can be accounted[fcr by the renormal- 
ization of the SO coupling constaniB3. Spin-orbit (SO) 
interaction emanates from the Pauli equation and is of 
the form: 



eh 



:V(p- (axp), 



(1.1) 



where me is the effective electron mass, (p is the scalar 
electromagnetic potential, a is the Pauli matrices vector, 
and p is the momentum operator. The vector electromag- 
netic potential A is put equal to zero in all the systems 
this paper deals with. 

The described system possesses a spectrum in which 
every energy level is twice degenerate. We introduce a 
perturbation (not necessarily small) to the confining po- 
tential that breaks axial symmetry of the system. This 
leads to lifting spin degeneracy, the mechanism of which 
being different depending on whether the perturbation 
has electromagnetic or non-electromagnetic origin. The 
investigation of effects, which such perturbations have 
on spectra, is yet more important if one takes into ac- 
count that perturbations conditioning deviations from 
axial symmetry are always present in the experiment. 
We show that under certain conditions the spin split- 
ting, originated by potentials in question, may be large 
enough to be detected. 



II. AXIALLY SYMMETRIC CONFINING 
POTENTIAL 



Let us first consider the case of axially symmetric con- 
fining potential created by purely electrostatic field. The 
inclusion of SO interaction into the Hamiltonian signifi- 
cantly hinders seeking the solution of the corresponding 
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eigenvalue problem, so that one may expect to find solu- 
tions only by means of perturbation theory or numerical 
calculations. That is why we chose a parabolic confining 
potential 



Vconf = , 



(2.1) 



where uj has the dimension of frequency and measures 
the strength of the potential and p is the lateral cylin- 
drical coordinate. This potential will provide us with 
a convenient zero approximation to the solution. The 
Hamiltonian of our problem has the following form: 



The radial functions ^ni (p) comprise the generalized La- 
guerre polynomials Ln . The value = ^/hjm^ serves 
as a characteristic scale in our problem. The function 
Rk^{z) corresponds to a plane wave with kz being the 
longitudinal wave number and normalization is based on 
the property S{0)/V ^ 1 as y ^ oo. The functions 
(2.6) form a complete orthonormal set. The energy cor- 
responding to the state Xnik^ is given by 



E, 



(0) 



huj{2n+ \l\ + 1) + 



21.2 



2mp 



(2.8) 



H 



2m, 



Here the SO interaction is 



2, ,2 



Vso 



iTl U! 



de 



2, ,2 



+ Vconf + PVSO + 



-ie- 
ie'" 



(2.2) 



d_ 

dz 



(2.3) 

The positive constant /i defines the strength of the SO in- 
teraction. At p — 1 the Hamiltonian (3) naturally stems 
from the usual Pauli Hamiltonian 



It is important that the whole energy may be separated 
into two conserving constituents: the lateral and longitu- 
dinal ones, and the lateral spectrum is equidistant. This 
property holds true as long as Hamiltonian depends on z 
only by the way of derivatives 3^ jdz^ . The longitudinal 
component of the wave vector, A:^, remains constant dur- 
ing the motion and is convenient to be considered as a 
parameter. At a given k^, the eigenfunctions (2.6) form a 
complete set in the corresponding subspace and may be 
chosen as a basis for seeking the approximate solutions 
of the eigenvalue problem 



Up = 



Imp 



- eif - 



P 



i7$ = E<^. 



(2.9) 



eh 



:a- {E xp) 



eh 



Smzc 



rdivE 



(2.4) 



if one puts the scalar potential ip equal to nif.uj'^p^ /2e and 
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neglects the term p /8m^c^, which does not influence the 
spin-splitting. 

The solution to the zero approximation eigenvalue 
problem 



1,2 

+ Vconf)x = Hox = Ex 



(2.5) 



was found in 1928 by FockEl. It reads 

XnlkAP,^,^) = ^ni{p)Yi{0)RkAz), (2.6) 



*ni(p) 



in+\l\)\pj'''^^ 2 (/J 



X 4(1 ( ) ,n = 0,1,2..., I = 0,±1,±2,..., (2.7) 



2tt 



The Hamiltonian (2.2) has a 2 x 2-matrix structure, so its 
basis eigenfunctions must be of a two-component vector 
form. Such two-component functions form the funda- 
mental system of the matrix Hamiltonian 



Ho 
Ho 



(2.10) 



and read: 



XT, nlk^ 



Q I ' Ai, nlk^ 





Xnlk^ 



(2.11) 



The diagonal (x(T,i), nikAVso\xa,i), mkkj and off-dia- 
gonal (X(t,x), „ifcjVso|X(i,T), mfcfcj matrix elements as 
well as matrix elements which appear in Section HI may 
be easily calculated. The convolutions of two-component 
vectors with 2x2 matrices were performed according 
to the usual matrix multiplication rules. The integrals 
needed to be taken in order to calculate the matrix ele- 
ments may be written as a product of three integrals over 
p, 9 and z. i^Ghe integrals over p were taken according to 
the formulae^: 
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{x)Xnl{x)dx 



r(l + a + + fc)/2) r(l + n + r((fc -l)/2-a + m) 



^m\n\ (TO+|fc|)!(n+|/|)! r(l + /) r(2 + (Z - fc)/2 + a) 



3F2 [-n, 2 + a + + /e)/2, 2 + a + (fc - 0/2; l + fc,2 + a- TO+(fc- 0/2; 1]. 



(2.12) 



In concordance with the normahzation rules for y;(6')E3 
the integrals over d and z are simple Fourier transfor- 
mations. For instance, what follows is a typical integral 
one meets in the calculation of the matrix elements from 
Section III: 



cos2(-l)e'('=-')«d0 = ^-{Uu + h-i,N + 5u-i.-n). 



'Y'^ 2 

(2.13) 

In numerical calculations throughout this paper we 
used 30 up- and 30 down-basis functions. This allowed 
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FIG. 1: (a) The spectrum of Ho; (b) the spectrum of H, 
/i = 1. Energy E is measured in eV, kz is measured in units 
of 10*^ cm-\ a = 8.47T0-^ eV. 



us to obtain 60 energy levels with a difference from the 
exact values only in the 9th significant digit. The value 
of Lo was chosen to be 10^^ sec^^ and m* was put equal 
to 0.05me. 

The Fig. ^ shows several first lateral energy levels of the 
Hamiltonians Hq and H. The lateral energies pertaining 
to Hq are degenerate with the multiplicity 2(2n-|- |Z| -I- 1). 
Thus, every g-th eigenvalue (counting from the bottom 
of the quantum well) of is 2q times degenerate and 
splits into q branches in transition to H . Therefore, Vso 
leaves 2-fold spin degeneracy of every energy level. 



III. PERTURBED AXIALLY SYMMETRIC 
CONFINING POTENTIAL 

A. Non-electromagnetic perturbation 

The perturbation we will deal with here and in the 
subsequent part of this section is 



Vpert — ^ 



A>0, iVG 0,1,2,, 



(3.1) 



(3.2) 



The dimensionless constant parameters A and N define 
the amplitude and the number of 'lobes' of the pertur- 
bation, respectively (see Fig. H). At = 2 our problem 
resembles the 2D type of the problem solved by MigdaS. 
In this part we investigate the effect which Vpert has on 
the spectrum considering it to be of non-electromagnetic 
nature {Vp~f ). Under the non-electromagnetic nature 
of the perturbation we imply that it does not affect Lp 
and E in the Hamiltonian (2.4), but is simply added to 
the right-hand side of (2.4). On the other hand, the elec- 
tromagnetic perturbation V^l^i , which is regarded in the 
next subsection, does influence ip and E. 

The case of the non-electromagnetic perturbation is es- 
pecially interesting because at even N V"~f leaves the 
inversion symmetry (0 — > -|- tt, z — > —z) of the Hamil- 
tonian intact, although breaks its axial symmetry. This 
should have a non-trivial effect on the spin splitting in 
the energy spectra. The increased (in comparison with 
the free space value n — 1) value of /z, equal to 10'^, ap- 
proximately corresponds to that of the constant /3 from 
the worm and may not only allow to draw correlations 
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FIG. 2: The polar plot of the potential Vconf + V^ert, A'^ = 6, 
p = puj'- dotted line - A = 1.0, dashed line - A — 0.5, solid 
line - A = 0.1. 



with the results of the latter, but can also help to eluci- 
date the role of SO interaction in the formation of energy 
band structure. The Fig. ^ demonstrates spin splitting 
for various values of the parameters N, A, and ^. At 
even values of N no spin splitting is predicted for any 
energy level. This may be explained by higher symmetry 
of the potential in this case. At odd N the spin splitting 
follows a rather peculiar pattern: with the growth of N 
the number of non-perturbed energy level (if the ground 
state has number zero) at which spin splitting occurs for 
the first time enhances. The Table 1 illustrates this be- 
havior. 



Table 1. Energy levels of the Hamiltonian H which 
undergo splitting in the case of non-electromagnetic per- 
turbation Vpert'^- The ground state has number zero. 



N 


1 


2 


3 


4 


5 


6 


7 


Level number 


all 


none 


2,3,9,... 


none 


4,6,... 


none 


7,10,... 



Very important for practical reasons is the question of 
how spin splitting depends on the amplitude of pertur- 
bation A. The Table 2 shows that the magnitude of spin 
splitting decreases very rapidly as A gets smaller, so that 
minute perturbations give almost indiscernible splitting. 
The rise of the contribution of SO interaction (the in- 
crease of fi) does not lead to any fundamental changes in 
the energy spectra, but increases almost proportionally 
to /i the values of spin splitting. For example, ai A = 1, 
iV = 3, = 103 and = 2-lO'^cm-i spin splittings \ Ad\, 
where d is the number of the unperturbed energy level, 
for d = 2 and d = 3 are 2.32 • 10"^ eV and 3.07 • 10"^ eV, 
respectively. 



1.820325 - 




..■■•**■ 


1.820320 - 






1.820315 - 






1.820310 - 






1.820305 - 




*"■•••.., 


1.820300 - 




"""■""^ 




' 1 ' 

5 


1 ■ 1 ' 1 ' 

10 15 20 



^1 

1.835 
1.830 
1.825 -I 
1.820 
1.815 
1.810 



(b) 



10 



15 



20 



(c) 



1.866 



1.864 



1.862 



1.860 



10 



15 



20 



FIG. 3: The spin splitting originated by non-electromagnetic 
perturbation V^"^'^'. Energy E is measured in eV, kz is mea- 
sured in units of lO^cm"^. (a) A = 1, /i = 1, = 3, d = 1, 2; 
(b) A = 1, p ^ 10^ iV = 3, d = 1,2; (c) yl = 1, /i = 10^ 
iV = 4, d= 1,2. 



B. Electromagnetic perturbation 



Treated as being of electromagnetic nature, the pertur- 
bation Vl'gVt' which formally has the form (3.1), results 
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in the following Hamiltonian: 

H' = H + AH. 



(3.3) 



Here 



AH - 



ArUpC^ 



Ap cos^( 




i)coii{Ne) 



AN ,Ne, . ,Ne, 

— pcos(— )sin(— ) 




AN ,N9 

pcosi 

2 ^ ^2 



. ,Ne, d 
sm(-— )(T, 



dp 



. 2, NO. d 



(3.4) 



r 



The matrix elements of AH in space of the func- 
tions (2.13) can also be explicitly calculated (see Eqns. 
(2.13)-(2.15)). The Fig. || demonstrates the spin splitting 
originated by V^eVt- Now all the energy levels are non- 
degenerate for both even and odd values of N . The mag- 
nitude of spin splitting goes down with the decrease of A 
approximately with the same rate as it does in the case 
of non-electromagnetic perturbation (see Table 2). It is 
clear from the Table 2 that as A decreases, the difference 
between non-electromagnetic and electromagnetic cases 
gradually vanishes. 



Table 2. The value of spin splitting (eV) at various 
amplitudes A of the perturbation Vpert- N = i, p ^ 1, 
k,=2 - W cm-i. 



A 


1 


10-1 


10-2 


10-3 


^rn — el 
^ pert 


IA2I 


2.29-10-^ 


4.16- 10-s 


4.06- lO-'^ 


4.05 - 10-8 


-Crn — el 
* pert 


IA3I 


3.04-10-5 


6.16- 10-<^ 


6.09-10-^ 


6.07-10-8 


* pert 


IA2I 


9.63-10-6 


2.87-10-6 


4.04-10-^ 


4.05-10-8 


* pert 


IA3I 


1.75-10-5 


4.65 - 10-6 


6.08-10-^ 


6.07-10-8 



IV. CONCLUSION 

In this work we studied theoretically a 3D gas of 
non-interacting electrons which was transformed into a 
QIDEG by means of 2D confining potential. We have 
solved the Pauli equation numerically for this system 
with perturbed axially-symmetric, parabolic confining 
potential. The non-perturbed axially symmetric poten- 
tial Vconf was assumed to be created by electromagnetic 
field with the scalar potential (p = meu'^ p'^ /2e. Two dif- 
ferent types of perturbations capable of lifting the spin 
degeneracy are considered: of electromagnetic and non- 
electromagnetic nature. We have found that these two 
types of perturbation have different effects on the energy 
spectra. While the electromagnetic perturbation splits 



all the levels, the non-electromagnetic one in most cases 
(besides the case with A'^ = 1) leaves the spin degeneracy 
of some levels. We have found that for large perturba- 
tions {A = 1) the value of spin splitting is large enough 
to be detected in the experiment (10-^ eV and higher). 
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FIG. 4: The spin splitting originated by electromagnetic per- 
turbation Vplrt- Energy E is measured in eV, kz is measured 
in units of 10® cm-^ (a) A = 1, = 1, iV = 3, d = 1, 2; (b) 



more complex behavior of high-lying levels, A — 1, p. 
iV = 4, d = 21, . . . , 25. 
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